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Abstract. One of the most striking results is the large elliptic flow (1)2) at RHIC. Detailed mass and 
transverse momentum dependence of elliptic flow are well described by ideal hydrodynamic calculations 
for pt < 1 GeV/c, and by parton coalescence/recombination picture for pt = 2 — 6 GeV/c. The systematic 
error on V2 is dominated by so-called "non-flow effects", which is the correlation not originated from 
reaction plane. It is crucial to understand and reduce the systematic error from non-flow effects in order 
to understand the underlying collision dynamics. In this paper, we present the centrality dependence of 
V2 with respect to the first harmonic event plane at ZDC-SMD (-^{ZDC-SMD}) in Au + Au collisions at 
y/SNN = 200 GeV. Large rapidity gap (\Arj\ > 6) between midrapidity and the ZDC could enable us to 
minimize possible non-flow contributions. We compare the results of H2{ZDC-SMD} with «2{BBC}, which 
is measured by event plane determined at \rj\ = 3.1 — 3.9. Possible non-flow contributions in those results 
will be discussed. 

PACS. 25.75.-.q Relativistic heavy-ion collisions - 25.75.Ld Collective flow 



1 Introduction 

Elliptic flow is expected to be one of the key observables 
to study an early stage of heavy ion collisions [T]. It is 
defined by the second harmonic Fourier coefficient 

W2 = <COB(2[0-!Pbp])>, (1) 

where (f> is the azimuthal angle of emitted particles, Wrf is 
the azimuthal angle of reaction plane and brackets denote 
the average over all particles and events. 

The PHENIX experiment at Relativistic Heavy Ion 
Collider (RHIC) have measured the v 2 for identified charged 
hadrons [2,3], <\> mesons and deuterons [4], 7r°'s and pho- 
tons [5] as well as electrons from heavy flavor decays [5] at 
midrapidity. The mass ordering of v 2 for identified hadrons 
were qualitatively explained by ideal hydrodynamics in 
transverse momentum < 2 GeV/c [2]. For intermediate 

= 2 — 6 GeV/c, a universal parton v 2 were obtained by 
dividing v 2 and px by constituent quarks for each hadron 
[2,3]. The v 2 for <p meson was also found to follow the 
quark number scaling, which support that the parton v 2 
have already developed prior to the hadronization [4 . Be- 
cause the cross section of (j) meson to non-strange hadrons 
are small, <j) meson v 2 is less sensitive to the late hadronic 
stage. The finite v 2 for electrons from heavy flavor decays 
implied the non-zero charm v 2 6] . Comparison of v 2 with 
transport model calculation suggest that the viscosity to 
entropy density ratio is close to the quantum lower bound 

1/47T 0. 



These measurements were done by using an event plane 
determined from the Beam-Beam Counter (BBC) located 
at pseudorapidity \r/\ — 3.1 - 3.9. The large pseudorapidity 
separation \Ar]\ ~ 3 from midrapidity would reduce non- 
flow effects, which are correlations not originated from the 
reaction plane such as jets, resonance decays and so on. 
Fluctuations of v 2 were also considered as the non-flow 
contributions 8], which would become more important in 
smaller system, such as Cu + Cu collisions. It is crucial 
to understand how non-flow contributions affect the event 
plane determined at the BBC so that we could validate a 
sensitivity to real collective flow on our measured v 2 based 
on the BBC. 

In this paper, we present the v 2 with respect to the 
event plane from directed flow determined at the Shower 
Maximum Detector (SMD), which is located at | | > 6. 
The larger rapidity separation could reduce the possible 
non-flow effects on our measured v 2 ■ We will compare the 
v 2 results from the event planes determined at the BBC 
and SMD and discuss the possible non-flow contributions 
on the v 2 . 



2 Data Analysis 

In this study, we analyzed ~ 650 M events collected by the 
PHENIX experiment in Au + Au at y/sWw = 200 GeV. 
Minimum bias events were selected within a collisions z- 
vertex ± 30 cm. Event centrality was determined by the 
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correlation between the energy deposit at the Zero Degree 
Calorimeter (ZDC) and number of charged particles at the 
BBC. Tracking were done by the Drift Chamber (DC) and 
Pad Chambers (PCs) at the central arm |ry| < 0.35. Trans- 
verse momentum were determined by the incident angle 
at the DC. The polar angle of the tracks were obtained by 
the hit at the inner PC (PCI) and the collision vertex from 
the BBC. Track associations were made by comparing hit 
positions with the projection of the DC tracks to the outer 
Pad Chamber (PC3). Tracks were required to have a hit 
on the PC3 within ± 3 a of the expected hit location in 
both azimuthal and beam directions. Large energy deposit 
E/p > 0.2 at the Electromagnetic Calorimeter (EMCal) 
were also required, where E is the energy deposit in EM- 
Cal and p is the momentum determined at the DC. It was 
necessary to eliminate the background mainly from pho- 
ton conversions and resonance decays so as to improve the 
signal to background ratio for pt > 4 GeV/c. 

The v 2 was measured by an event plane method [5] and 
was obtained by dividing the measured v 2 by the event 
plane resolution 



,,obs 



C 2 



(cos (2[</>-V n ])) 

Res{^ n } (cos(2[tf n -% P ]))- 



(2) 



where <fi is the azimuth of charged hadrons at the central 
arm (|?7| < 0.35), !?„ is the event plane from the n-th 
harmonic flow (n = 1 for the ZDC-SMD, n = 2 for the 
BBC) and v 2 bs is the measured v 2 with respect to the event 
plane tfv Event planes were determined from the v 2 at the 
BBC and the central arm as well as the directed flow v\ — 
(cos (<j) — ■Z'rp)) at the Shower Maximum Detector (SMD). 
The central arm event plane is only used to evaluate the 
event plane resolutions. The SMDs are located at the same 
acceptance of the ZDCs, \rj\ > 6, and measure transverse 
positions of spectator neutrons. The measured v 2 s are 
denoted as w 2 {ZDC-SMD} and v 2 {BBC} for the ZDC- 
SMD and BBC event planes, respectively. 

The event plane determined at the ZDC-SMD can min- 
imize non-flow correlations as well as v 2 fluctuations be- 
cause of the following reasons. First, the pseudorapidity 
gap from midrapidty is 6, which is higher than what we 
have previously studied the v 2 by using the BBC. Second, 
the ZDC-SMD event plane is determined from directed 
flow. This mixed harmonic method involves three parti- 
cle correlations and thus direct two particle correlations, 
which is dominant contributions from non-flow effects, do 
not affect the measured v 2 . Third, ZDC-SMD measures 
spectator neutrons rather than participants. Therefore, v 2 
fluctuations are suppressed up to the fluctuation of spec- 
tator neutrons. 

Fig. [1] shows the event plane resolutions as a function 
of centrality. At least two independent event planes are 
required in order to evaluate the resolution since the az- 
imuth of true reaction plane is unknown. The resolution 
from two independent event planes is calculated by 



Res{^} = C J {cos (2[Wn - #+])), 



(3) 
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Fig. 1. Event plane resolutions as a function of centrality for 
the ZDC-SMD event plane (solid circles) and the BBC (open 
diamonds) by Eq. |3j. Dashed lines represent the resolutions 
calculated by Eq. (f4]| . 



A constant parameter C is very close to V2 for both BBC 
and ZDC-SMD due to low resolution [9]. The ZDC-SMD 
resolution is about factor 4 smaller than that of BBC be- 
cause the ZDC-SMD event plane is determined from di- 
rected flow. It is approximately proportional to u^Af SMD , 
where M SMD is multiplicity used to determine the ZDC- 
SMD event plane, whereas the BBC resolution is roughly 
proportional to v 2 V M BBC . 

The resolutions were also evaluated by adding a refer- 
ence event plane 



Res{#„} = 



/ (co S (2[^-^ n ]))(cos(2[^ n -^])) 
<cos(2[<^-^])) 



where and \P n denote the event planes determined at 
the forward and backward pseudorapidities, respectively. 



( 4 ) 

where I, m and n denote the harmonics for event plane 
ty A , 'P B and respectively. Dashed lines in Fig. [1] show 
the resolutions calculated by Eq. (T4J. For example, the 
BBC resolution was calculated by inserting fy n — ty® 30 , 

W A = ^ZDC-SMD and = ^CNT where CNT denote 

the central arm. One can find that the dashed lines are 
systematically lower for the ZDC-SMD, and higher for the 
BBC. The comparison of v 2 from two different resolutions 
will be presented in the next section. 



3 Results 



We will present the preliminary results of ^{BBC} as 
well as {ZDC-SMD} in Au + Au collisions at ^/sjvtv = 
200 GeV measured at the PHENIX experiment. Section 
13. II will give comparison of the v 2 between BBC and ZDC- 
SMD event planes. Results between PHENIX and STAR 
experiments will be compared in Section 13.21 In Section 
13.31 centrality dependence of the w 2 {ZDC-SMD} will be 
compared with the D2{BBC}. 
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Fig. 2. (Top) Comparison of V2 as a function of pt in 20 - 
60 % centrality for the BBC (open diamonds), the ZDC-SMD 
from two different event plane resolutions (open crosses and 
solid circles, see texts). Gray bands, solid red lines and yel- 
low boxes represent systematic uncertainties on the U2{BBC} 
and v 2 {ZDC-SMD}. (Bottom) The ratio of u 2 {ZDC-SMD} to 
«2{BBC} as a function of px- Dashed line denote the fit result 
by constant. 



3.1 Comparison of v 2 {BBC} with v 2 {ZDC-SMD} 

Fig.dlshows the u 2 {ZDC-SMD} as a function of p T in 20 
- 60 % centrality. For comparison, the u 2 {BBC} is also 
plotted by open diamonds. The v 2 increases linearly up to 
Pt ~ 3 GeV/c, reach maximum ^0.2 and then start de- 
creasing for higher p T . The v 2 {ZDC-SMD} (S-N), which 
is obtained from the resolution in Eq. ((3]), is about 7 % 
systematically lower than the w 2 {BBC}, while the results 
are consistent within systematic uncertainties. We also 
plot the w 2 {ZDC-SMD} (ZDC-BBCS-BBCN) as shown 
by solid circles, which is obtained from the resolution in 
Eq. (0} by inserting <P n = ^ ZDC - SMD , = <Pf BCS and 
&z = ^bbcn The BBCg and BBCN denote the backward 

and forward BBC, respectively. The w 2 {ZDC-SMD} from 
two different resolutions are in good agreement within sys- 
tematic uncertainties. Bottom panel shows the ratio of 
u 2 {ZDC-SMD} to w 2 {BBC} as a function of p T . One can 
see that the ratio is constant within systematic errors in 
the measured pt range. 



3.2 PHENIX vs STAR 

Fig. [3] show the comparison of the PHENIX ^{ZDC- 
SMD} with STAR result [10] in 20 - 60 % centrality bin. 
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Fig. 3. Comparison of the u 2 {ZDC-SMD} between PHENIX 
and STAR experiments in 20 - 60 % centrality. The STAR 
«2{ZDC-SMD} is taken from |10| . Gray bands and light blue 
boxes denote the systematic uncertainties on the «2{BBC} and 
«2{ZDC-SMD}, respectively. Bottom panel shows the ratio of 
V2 to the fourth polynomial fit of the «2{BBC} as a function 
of pt. 



Only statistical errors are shown for the STAR t> 2 . Both 
PHENIX and STAR results are obtained by the resolution 
in Eq. © and thus the results of ^{ZDC-SMD} are ex- 
tracted by the exactly same method. Data symbols (open 
diamonds and open crosses) are the same as shown in 
Fig. O For a quantitative comparison, the ratio of v 2 to 
the f 2 {BBC} is plotted in bottom panel in Fig. [3] The de- 
nominator of the ratio is fitting result of the w 2 {BBC} by 
fourth polynomial function. One can see that the results 
agree very well within systematic errors. 



3.3 Centrality dependence of w 2 

Fig. 0] show the v 2 {pt) for 10 % step centrality bin in 20 
- 60 %. In left figures, the u 2 {ZDC-SMD} is essentially 
consistent with the 1> 2 {BBC} within systematic errors. In 
peripheral 40 - 60 %, we find that the u 2 {ZDC-SMD} is 
5 - 10 % lower than the v 2 {BBC}. The lower ^{ZDC- 
SMD} could suggest possible non-flow contributions on 
the w 2 {BBC}. 

In right figures, both BBC and ZDC-SMD resolutions 
are calculated by ZDC-SMD, BBC and CNT combina- 
tions in Eq. ([4}. Since the non-flow effects are expected 
to be maximum at the midrapidity, the CNT event plane 
could have maximal sensitivity to the non-flow contribu- 
tions. Therefore, by including the CNT event plane reso- 
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Fig. 4. Comparison of «2{ZDC-SMD} to ^{BBC} as a function of pt for 10 % step centrality bin in 20 - 60 %. Left (a) and 
right (b) figures show the results without and with the central arm (CNT) event plane resolution, respectively. 



lution, we could study how non-flow contributions on the 
CNT resolution modify the v 2 {BBC} as well as ^{ZDC- 
SMD}. We find that the v 2 {ZDC-SMD} become more con- 
sistent with the «2{BBC} by including the CNT event 
plane resolution. Because the ZDC-SMD (BBC) resolution 
from the ZDC-BBC-CNT combination in Eq. (@J is lower 
(higher) than that from Eq. ((3]) as shown in Fig. [T] Re- 
sulting w 2 {ZDC-SMD} (v 2 {BBC}) become higher (lower) 
than those from the resolution by forward and backward 
correlations. Therefore, the ^{ZDC-SMD} is closer to the 
U2{BBC} when one include CNT event plane to evaluate 
the resolution. 



4 Conclusion 

In summary, we have measured unidentified charged hadron 
elliptic flow with respect to the ZDC-SMD event plane 
from directed flow in Au + Au at ^/snn = 200 GeV. 
The i>2{ZDC-SMD} was compared with the v 2 measured 
with respect to the event plane determined at the BBC. 
We found that the ^{ZDC-SMD} was basically consis- 
tent with the «2{BBC} within systematic uncertainties in 
20 - 60 % centrality. Several different choice of event plane 
resolutions were studied. We found that resulting w 2 {ZDC- 
SMD} was still consistent with the «2{BBC} even if the 
CNT event plane was included in the event plane resolu- 
tion. The difference of v 2 between BBC and ZDC-SMD 
event plane, ~ 5 - 10 %, at 40 - 60 % centrality bins could 
attribute to the possible non-flow effects. This result in- 
dicate that the non-flow effects are essentially minimal on 
the v 2 {BBC} in 20 - 60 % centrality because the ^{ZDC- 
SMD} is expected to be unbiased by the non-flow contri- 
butions. 
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